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Monitoring the fluorescent radiation of an atom unravels the master equation evolution by col- 
lapsing the atomic state into a pure state which evolves stochastically. A robust unraveling is one 
that gives pure states that, on average, are relatively unaffected by the master equation evolution 
(which applies once the monitoring ceases). The ensemble of pure states arising from the maximally 
robust unraveling has been suggested to be the most natural way of representing the system [H.M. 
Wiseman and J. A. Vaccaro, Phys. Lett. A 250, 241 (1998)]. We find that the maximally robust un- 
raveling of a resonantly driven atom requires an adaptive interferometric measurement proposed by 
Wiseman and Toombes [Phys. Rev. A 60, 2474 (1999)]. The resultant ensemble consists of just two 
pure states which, in the high driving limit, are close to the eigenstates of the driving Hamiltonian 
Qa x /2. This ensemble is the closest thing to a classical limit for a strongly driven atom. We also 
find that it is possible to reasonably approximate this ensemble using just homodyne detection, an 
example of a continuous Markovian unraveling. This has implications for other systems, for which 
it may be necessary in practice to consider only continuous Markovian unravelings. 
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I. INTRODUCTION 

Some states of open quantum systems are more ro- 
bust than others. That is, they are less perturbed by 
the system dynamics. This fact has been the subject of 
a long-running and active research program pj-Q]. Re- 
cently, one of us and Vaccaro have introduced into this 
program a formalism with a number of distinctive fea- 
tures |^,^]. This formalism consists of finding the max- 
imally robust unraveling (MRU) for the open quantum 
system. Its introduction was motivated by a desire to 
better understand the rich dynamics of open quantum 
systems in general || , and that of the atom laser in par- 
ticular |J. 

The use of the term "robust unraveling" rather than 
"robust state" encapsulates two of the distinctive features 
of the work of Refs. . An unraveling is a way of mea- 
suring the environment of an open quantum system such 
that the system state can be described by a pure state 
undergoing stochastic evolution. This is always possible 
in principle if the unmonitored system obeys a Marko- 
vian master equation, which we will assume to have a 
unique stationary state. In the long time limit, the "un- 
raveled" system will be in a pure state drawn at random 
from a particular ensemble of pure states defined by the 
unraveling. It is from the consideration of such an en- 
semble of pure states that the two distinctive features of 
our approach are met. The first is that it is not indi- 
vidual pure states whose robustness are to be calculated, 
but rather a whole ensemble of pure states, the average 
robustness of which is calculated. The second is that the 
pure states in this ensemble are physically realizable in 
the sense that they are the states of the system known 
to an experimenter using the appropriate measurement 
scheme on the system's environment. 

As noted above, the principle application of the max- 
imally robust unraveling formalism has been to a model 



for an atom laser (a continuously damped and replen- 
ished gaseous Bose-Einstein condensate). This work is 
a specialized application in two ways. First, the system 
itself has a high excitation number and so is a quantum 
system in the classical limit. Second, only a subset of 
the set of all possible unravelings was considered. This 
subset (which is still infinite) contains those unravelings 
that lead to continuous and Markovian evolution of the 
system state vector [ fLOj . This restriction was necessary 
to make the problem tractable and was justified by the 
classicality of the system. 

In this work we apply the formalism of MRU to a sys- 
tem with no classical limit (in the usual sense at least), 
a resonantly-driven fluorescent two-level atom. This sys- 
tem is one of the canonical examples of an open quan- 
tum system, and has surprisingly complex dynamics for 
its size. It is therefore worth investigating in its own 
right. But, even more importantly, it is simple enough 
that the maximally robust unraveling can be found ana- 
lytically. It turns out that this MRU is neither continuous 
nor Markovian. This enables us to investigate the ques- 
tion of how closely one can approximate the ensemble of 
this MRU if one is restricted to considering continuous 
Markovian unravelings. The answer to this question has 
implications for the general usefulness of the MRU for- 
malism, since for typical systems it would be necessary 
to impose this restriction in order to make the formalism 
practical. 

The structure of this paper is as follows. In Sec. II we 
briefly review the MRU formalism. In Sec. Ill we intro- 
duce the two-level atom model and derive an expression 
for the ensemble average survival probability (which is 
used to quantify the robustness) in terms of moments of 
the ensemble of state vectors. In Sec. IV we look at one 
simple (but not maximally robust) unraveling, that re- 
sulting from direct detection of the atom's fluorescence, 
for comparison with other, more robust, unravelings. In 
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Sec. V we present the most robust unraveling and its en- 
semble of (in this case, just two) state vectors. In Sec. VI 
we find the most robust unraveling from within the set 
of continuous Markovian unravelings. We compare this 
ensemble to the MRU of Sec. V in Sec. VII. We con- 
clude with a discussion of the implications of our results 
in Sec. VIII. 



II. MAXIMALLY ROBUST UNRAVELINGS 
A. The Master Equation 



where the \ipk) are normalized state vectors and the pk 
are positive weights summing to unity. The (possibly 
infinite) set of ordered pairs, 

E = {(\1> k ),p k ):k}, (2.5) 

we will call an ensemble E of pure states. Note that there 
is no restriction that the states be mutually orthogonal. 
This means that there are continuously infinitely many 
ensembles E that represent p ss . The aim of finding the 
MRU is to find the "best" or "most natural" representa- 
tion for p ss . 



Open quantum systems generally become entangled 
with their environment, and this causes their state to 
become mixed. In many cases, the system will reach an 
equilibrium mixed state in the long time limit. This is 
the sort of system for which our approach to robustness, 
of finding the maximally robust unraveling (MRU), can 
be applied without modification. 

If the system is weakly coupled to the environmental 
reservoir, and many modes of the reservoir are roughly 
equally affected by the system, then one can make the 
Born and Markov approximations in describing the ef- 
fect of the environment on the system Jll| . Tracing over 
(that is, ignoring) the state of the environment leads to 
a Markovian evolution equation for the state matrix p of 
the system, known as a quantum master equation. The 
most general form of the quantum master equation that 
is mathematically valid is the Lindblad form (l^] 

M 

p=-i[H,p} + Y / Vlc f t}p = £p, (2.1) 

where for arbitrary operators A and B, 

V[A]B = ABA* - (A^AB + BA^A) /2. (2.2) 

If the master equation has a unique stationary state 
(as we will assume it does), then that is defined by 

Cp ss = 0. (2.3) 

This assumption requires that L be time- independent. In 
many quantum optical situations, such as resonance flu- 
orescence, one is only interested in the dynamics in the 
interaction picture, in which the free evolution at optical 
frequencies is removed from the state matrix. Indeed, if 
one treats the driving field as classical, as we will do, it 
is necessary to move into such an interaction picture in 
order to obtain a time-independent Liouvillian superop- 
erator C. 

The stationary state matrix p ss can be expressed as an 
ensemble of pure states as follows: 

p ss = Y,Pk\ipk)(4>k\, (2.4) 

fc 



B. Unravelings 

As explained in the Introduction above, the first crite- 
rion for our most natural ensemble is that it be physically 
realizable by monitoring the environment of the system. 
In the situation where a Markovian master equation can 
be derived, it is possible (in principle) to continually mea- 
sure the state of the environment on a time scale large 
compared to the reservoir correlation time but small com- 
pared to the response time of the system. This effectively 
continuous measurement is what we mean by "monitor- 
ing" . In such systems, monitoring the environment does 
not disrupt the system-reservoir coupling and the system 
will continue to evolve according to the master equation 
if one ignores the results of the monitoring. 

By contrast, if one does take note of the results of mon- 
itoring the environment, then the system will no longer 
obey the master equation. Because the system-reservoir 
coupling causes the reservoir to become entangled with 
the system, measuring the former's state produces infor- 
mation about the latter's state. This will tend to undo 
the increase in the mixedness of the system's state caused 
by the coupling. 

If one is able to make perfect rank-one projective (i.e. 
von Neumann) measurements of the reservoir state, the 
system state will usually be collapsed towards a pure 
state. However this is not a process that itself can be 
described by projective measurements on the system, be- 
cause the system is not being directly measured. Rather, 
the monitoring of the environment leads to a gradual (on 
average) decrease in the system's entropy. 

If the system is initially in a pure state then, under 
perfect monitoring of its environment, it will remain in a 
pure state. Then the effect of the monitoring is to cause 
the system to change its pure state in a stochastic and (in 
general) nonlinear way. Such evolution has been called 
a quantum trajectory p3| , and can be described by a 
nonlinear stochastic Schrodinger equation for the system 
state vector Q-Q. The nonlinearity and stochasticity 
are present because they are a fundamental part of mea- 
surement in quantum mechanics. 

On average, the system still obeys the master equa- 
tion. That is, if the increment in \if;} under the SSE is 
\di/j) then 
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E[|#)(V| + + - £E[|V»>W (2-6) 



^W = (^|e £t [|^><V#I] IV#> 



(2.9) 



Here E denotes the ensemble average with respect to the 
stochasticity of the SSE. This stochasticity is evidenced 
by the necessity of retaining the ltd term \dip)(dip\ fl7f| . 

Because the ensemble average of the system still obeys 
the master equation, the stochastic Schrodinger equa- 
tion) is said to unravel the master equation [jl3|. It is 
now well-known that there are many (in fact continuously 
many) different unravelings for a given master equation 
Ji8|] , corresponding to different ways of monitoring the 
environment. 

Each unraveling 11 gives rise to an ensemble of pure 
states 



{(h#>, *#):*}, 



(2.7) 



where are the possible pure states of the system 

at steady state, and p^ are their weights. For master 
equations with a unique stationary state p ss , the SSE is 
ergodic over E u [jl9| and p" is equal to the proportion 
of time the system spends in state The ensemble 

E u represents p ss in that 



k 



(2.8) 



as guaranteed by Eq. (|2.6| ). 



C. Survival Probability 

Imagine that the system has been evolving under a 
particular unraveling 11 from an initial state at time — oo 
to the stationary ensemble at the present time 0. It will 
then be in the state \tp^) with probability p^ . If we now 
cease to monitor the system then the state will no longer 
remain pure, but rather will relax toward p ss under the 



evolution of Eq. ( |2.f| ) . 

This relaxation to equilibrium will occur at different 
rates for different states. For example, some unravel- 
ings will tend to collapse the system into a pure state 
that is very fragile, in that it changes into a very dif- 
ferent (and mixed) state as it relaxes to equilibrium. In 
this case the ensemble would rapidly become a poor rep- 
resentation of the observer's expected knowledge about 
the system. Hence we can say that such an ensemble is 
a "bad" or "unnatural" representation of p. Conversely, 
an unraveling that produces robust states would remain 
an accurate description for a relatively long time. We ex- 
pect such a "good" or "natural" ensemble to give more 
intuition about the dynamics of the system. The most 
robust ensemble we interpret as the "best" or "most nat- 
ural" such ensemble. 

We quantify the robustness of a particular state 
by its survival probability S^(t). This is the probability 
that the system would be found (by a hypothetical pro- 
jective measurement) to still be in the state at time 
t. It is given by 



Since we are considering an ensemble E u we must de- 
fine the average survival probability 



(2.10) 



In the limit t — > oo the ensemble-averaged survival prob- 
ability will tend towards the stationary value 



S w (oo)=Tr[p 8 2 s ]. 



(2.11) 



This is independent of the unraveling 11 and is a measure 
of the mixedness of p ss . 

There are many possible figures-of-merit that may be 
obtained from the survival probability S u (t), as discussed 
in Rcf. H . Here we choose the simplest one, also adopted 
in P): the time it takes for S u (t) to fall half-way to its 
equilibrium value. That is, 



min{i : S U (t) = (1 + Tr[p 2 s ])/2}. 



(2.12) 



This survival time r w quantifies the robustness of a par- 
ticular unraveling 11. 

Let the set of all unravelings be denoted J. Then the 
subset of maximally robust unravelings Jm is 



J M = {1Z E J : r n > t u VU e J}. 



(2.13) 



Even if Jm has many elements 72-1,72.2, ■ ■ ■> these different 
unravelings may give the same ensemble E n = E ni = 
E U2 = . . .. In this case we claim E n is the most natu- 
ral ensemble representation of the stationary solution of 
a given master equation. Different definitions of survival 
time pH will obviously lead to different numerical values 
for t . We are less concerned with such numerical val- 
ues than with the robust ensemble E 1 ^ , which has been 
found p| to depend little on the precise definition used. 



III. THE TWO-LEVEL ATOM 

A. The Resonance Fluorescence Master Equation 

Consider an atom with two relevant levels {|g),|e)}. 
Let there be a dipole moment between these levels so 
that the coupling to the continuum of electromagnetic 
field modes in the vacuum state will cause the atom to 
decay at rate 7. So that the atom does not simply decay 
to the state \g), add driving by a classical field (such as 
that produced by a laser) of Rabi frequency ft. We work 
in the interaction picture with respect to the free Hamil- 
tonian Hq — huio\e)(e\ so that the classical driving at 
frequency to becomes time-independent. The evolution 
of the atom's state matrix can then be described by the 
resonance fluorescence (RF) master equation 
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P = -i^[<7x,p\ +lV Hp. 



(3.1) 



In this equation we have used the Pauli matrices 



(3.2) 
(3.3) 
(3.4) 

(3.5) 
(3.6) 



In terms of these, any state of the atom can be written 
as a 3-vector (x,y,z) satisfying 



<Jx = 


\e)(g\ + \g){e\ 


(Ty = 


-i\e)(g\ + i\g)( 


a z = 


\e)(e\-\g)(g\ 






a = 


\g)(e\ = - 


at = 


\e)(g\ = 



x 2 + y 2 + z 2 < 1, 



(3.7) 



with equality only for a pure state. From this Bloch vec- 
tor the state matrix is defined by 



1 , 

p= - (I + xa x + yo-y + za z ) . 



(3.8) 



The linear equatio ns o f motion for the Bloch vector 
that result from Eq. (3.1) are known as the Bloch equa- 
tions. The solution satisfying the initial condition 



x(0) = u , y(0) = v , z(0) = w 



x(t) = ue-(~<W t , 

y(t) = c+e A+t +c_e A -' 



Vs, 



z(t) = c+ 



An m 

Here c± are constants given by 
1 



+ Zm 



8in 



=F4fi(w - z ss ) ± (7 ± 4ifi)(t> - y s: 



The eigenvalues A± are defined by 



X± = ~^7 i *fi< 



Here 



n = v^ 2 - (7/4) 2 



(3.9) 



(3.10) 
(3.11) 

(3.12) 



(3.13) 



(3.14) 



(3.15) 



is a real modified Rabi frequency for £1 > 7/4, and is 
imaginary for f2 < 7/4. The stationary solutions appear- 
ing in the above equations are 



Us 



0. 



2 7 



7 2 + 29? ' 

~7 2 
7 2 + 2fi 2 ' 



(3.16) 
(3.17) 

(3.18) 



B. The Survival Probability 

Using the Bloch vector representation of the atomic 
state matrix it is easy to show that the survival proba- 
bility for a pure state \ip k ) with projector 
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\4>k){ipk\ = 7: (1 + U k<?x + V k CTy + w k a z ) 



(3.19) 



Sk(t) = - (1 + x k (t)u k + y k (t)v k + z k (t)w k ) , (3.20) 

where (x,y, z) k (t) is the Bloch vector at time t with 
the i nitial condi tion (x,y,z) k (0) — (u,v,w) k as in 
Eqs. ( 3.1Cl| )— ( 3.12 ). From that solution it is evident that 
S k (t) will contain terms that are constant, linear, and 
bilinear in the vector components (u, v, w) k of the initial 
state. 

As explained in the preceding section we are interested 
in the survival probability not for a single state but for 
an e nsem ble of states. This is the ensemble average of 
Eq. j2~10| ). After some work, the survival probability in 
this case is found to be simply 



S{t) = $>kS*(t) 



(3.21) 



i(l + y s 2 s + 4) + i 



(3.22) 

(1 - y 2 ss - 4)e-^/ 2 )' + V v f+(t) + V m f-(t) 



where 



/ ± ( t ) = _ e -(7/2)t + e -(37/W 



cos Qt ± — =- sin fit 



(3.23) 



In Eq. ( 3.22 ), all the information about the ensemble is 
contained in the moments 



V v = E[v 2 } - E[v} 2 = J2 PA - (j2 Pkv k ) , 

k \ k ) 

V w = E[w 2 } - E[w} 2 = PkW 2 k - ^ kWk I 

k \ k / 



(3.24) 



(3.25) 



This is possible because we have used the following rela- 
tions: 



E[(u,v,w)] = (x ss ,y ss ,z ss ), 

E[w 2 ] = 1 -E[v 2 ] - E[w 2 ]. 



(3.26) 
(3.27) 



To find the robustness of any particular unraveling we 
thus need to find simply the two ensemble averages V v 
and Vv,. 
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IV. UNRAVELING BY DIRECT DETECTION 

The most obvious way to unravel the RF master equa- 
tion is by direct detection. This requires detecting all of 
the atom's fluorescence by unit-efficiency photodetectors. 
This is beyond current technology, but not by so much 
that the experiment should be considered unphysical. As 
we will find, unraveling by direct detection is actually not 
very robust (by the definition of Sec. II) but it is never- 
theless useful to consider as a point of comparison with 
more robust unravelings. 

The stochastic evolution of an atom undergoing RF 
with direct detection has been considered many times be- 
fore . It has one feature that enables an enormous 
simplification over a generic unraveling. This is that im- 
mediately following a detection, the atomic state is inde- 
pendent of its state before, and is just the ground state 
\g). Between these jumps to the ground state, the condi- 
tioned atomic state evolves deterministically. At steady 
state, when there has certainly been at least one detec- 
tion, all members of the ensemble are therefore identified 
simply by the time t since the last detection. 

If there is a detection at time to, then, until the next 
detection occurs, the state of the atom at time to + t 
evolves according to the equation 13 R]] 



|l*,(*)> 



7 t 
—a' a — i — (j i 

2 2 



\Mt))- 



(4.1) 



The solution, satisfying the initial condition |-0o(O)) = \g) 
is 



where 

c„(*) 

Ce(t) 

Here 



|^o> =c e (t)\e)+c g (t)\g), 



cos {tit/ 2) + ^sm{Clt/2) 
-^sin(^/2)e- (7/4) *. 



,-(7/4)i 



n = v 7 ^ 2 - (7/2) 5 



(4.2) 

(4.3) 
(4.4) 

(4.5) 



is a real modified Rabi frequency for Q > 7/2 and is 
imaginary for £ 1 < 7 /2. Note that it is different from (l 
defined in Eq. (fj|). 

The state in Eq. (^^) is unnormalized, and the norm 
(V'o(t)lV'o(t)) represents the probability that there has 
been no detection since time to, given that there was 
a detection at that time. Let us write this probability as 



P {t) = \c e {t)\ 2 + \c g {t)f 



(4.6) 



We show in the appendix that this probability is related 
to pit), the probability that, at steady state, the last 
detection was a time t ago, by 



Pit) 



Po(t) 



f °°P {s)ds 



(4.7) 



As noted above, in steady state under direct detection 
the possible atomic states are parametrized by the real 
variable t, the time since the last detection. The state at 
that time has projector 



Pit) 



\Mt))(Mt)\ 
(Mt)\Mt)) ' 



(4.8) 



and the weight for each of these members of the ensemble 
is p{t)dt. Physically, all members of the ensemble exist 
on the u = great circle of the Bloch sphere , be cause that 
is where the modified Rabi cycling of Eq. ( f4.l| ) takes the 
ground state. This distribution is shown in Fig. [l](a). 
For f2 imaginary (that is, £1 < 7/2) the states never 
reach the excited state. For (l real (that is, £1 > 7/2, 
the states may undergo an arbitrary number of cycles. 
For !1>7 the states are likely to undergo many cycles 
before a spontaneous emission event occurs so that the 
ensemble consists of all the states on the u = great 
circle, almost un iformly distributed. 
From Eq. (L_ 



it can be verified analytically that 



P{t)p{t)dt = 



dt\Mt))(Mt)\ = p ss , (4.9) 



where 



[\~c e {t)f + \c a {t)\ 2 ] dt- 



(4.10) 



. Moreover we can easily find numerically the ensemble 
averages necessary to find the ensemble average survival 
probability, namely 



V v 



1 

a7 



[iC e (t)c* g (t) - iC* e (t)Cg(t)]' 

\c e {tW + \c g {t)\ 2 



dt 



S)f 



\c e {tW + \c g {tW 



-dt 



(4.11) 
(4.12) 



In Fig. g we plot the survival probability for the direct 
detection ensemble, for a variety of driving strengths Q. 
We see that for 11 < 7 the survival probability decays ap- 
proximately exponentially at rate of order 7. For fi -C 7 
the stationary state matrix is close to the ground state, 
and most members of the direct detection ensemble are 
also. For fi 3> 7 the ensemble is equally spread over 
the u = great circle and consequently the variances V v 
and V w are approximately equal to 1/2. Using this, the 
survival probability is found to be approximately 

S{t) « ~ (l + e -( 3 / 4 )T* C os tit) . (4.13) 

The oscillations in the survival probability are due to 
the Rabi oscillations. As noted above, the direct detec- 
tion ensemble consists of states on the u = great cir- 
cle. Rabi cycling around the a;-axis according to the RF 
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master equation ( p.l[ ) rotates this circle around, rapidly 
moving the states away from their initial positions and 
then back close to their initial conditions after one cycle. 
They do not return exactly to their initial states because 
of the slow (at rate 37 /4) decay towards the equilibrium 
state. This behaviour is illustrated for a typical member 
of the direct detection ensemble in Fig. |l|(b). 

The change from damped to oscillatory behaviour has 
a dramatic effect on the survival time in Eq. ( p. 12 ). It 
is plotted in Fig. || as a function of il. For £1 <§C 7 it is 
given by 



2 In 2 7 - 



(4.14) 



as in this limit the survival prob ability decays as e~ 7 */ 2 . 
For fl 7, we can use Eq. ( [4.13 ) to get the approximate 
expression 



-Q- 



(4.15) 



That is, the survival time is here determined by the 
Hamiltonian evolution only. As shown in Fig. |3j, this 
is quite a good approximation even for moderate fl. 



V. THE MOST ROBUST UNRAVELING 
A. The Most Robust Ensemble 

It was stated above that the unraveling by direct de- 
tection is not the most robust unraveling. In fact, from 
an examination of the survival probability in Eq. ( 3.22| ) 
we can see that it is one of the least robust unravclings 
for 3> 7. That is because of the large variances in v 
and w in this limit. 

It is not difficult to sho w tha t the two functions f+(t) 
and /-(t), defined in Eq. ( [3.23] ), are non-positive for all 
f2 and for t > 0. Since the variances V v and V w are non- 
negative, it is easy to see that to maximize the survival 
probability, one would wish to minimize V v and V w . The 
ideal limit would be V v = V w = 0. This corresponds to 
an ensemble in which all members have the same Bloch 
vector components v and w. Since the ensemble average 
must equal p ss it follows then that for all members 



(5.1) 



Furthermore, since x ss = 0, and since the members of the 
ensemble must be pure, it follows that for all members 
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(5.2) 



where the two alternatives are equally weighted. These 
two members of the ensemble are shown in Fig. 0(c). 

This ensemble is guaranteed to give the maximum sur- 
vival probability 



S(t) = -(l + y 2 ss + z 2 ss ) + - (l-i/, 



(5.3) 



which is plotted in Fig. |^. It gives the maximum survival 
time 



T = 2L127- 1 , 



(5.4) 



which is independent of Q. There is no Rabi cycling be- 
cause the states defined by the Bloch vector (u, v, w) have 
v and w already equal to their stationary values. Under 
the master equation evolution u simply decays towards 
its stationary value of zero and v and w remain constant. 
This decay of the Bloch vector to equilibrium at rate 7/2 
is shown in Fig. |l|( 



B. Adaptive Interferometric Detection 



The most robust ensemble defined here would be a 
mere curiosity if it were not for the fact that there is a 
detection scheme that realizes it. This scheme, proposed 
by one of us and Toombes [ po| , involves interfering the 
light from the atom with a resonant local oscillator before 
detection. This is done using a highly transmitting beam 
splitter as shown in Fig. ^. In the limit of a large local os- 
cillator, this is known as homodyne detection. However 
we require a very weak local oscillator, with reflected in- 
tensity comparable to the intensity of the light from the 
atom. Furthermore, we require the local oscillator ampli- 
tude to be continually adjusted by a real-time feedback 
loop. To be specific, the field detected should be propor- 
tional to 



(5.5) 



where the field from the atom is proportional to a, the 
atomic lowering operator as usual, and the local oscilla- 
tor field is represented by the complex number fj.(t). This 
complex number is given by 



/*(*) = ±: 



1 



(5.6) 



where the sign is changed every time a detection occurs. 

Remarkably, this relatively simple detection scheme 
has the consequence that, after initial transients have 
died away, a driven atom jumps between the states with 
projector P±, where 



Zs S (J x + v ss ay + w ss a z ) , (5.7) 



every time a detection occurs. The rate of detections 
moreover is independent of which of these states the atom 
is in, and is equal to 7/4. Thus in the long time limit the 
atom will have a probability p± — 1/2 to be in each of 
them, and the maximally robust ensemble will be physi- 
cally realized. 
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VI. CONTINUOUS MARKOVIAN 
UNRAVELINGS 

The most robust unraveling is neither continuous nor 
Markovian. It is not continuous because the atomic state 
jumps every time there is a detection. It is not Marko- 
vian because the evolution of the atomic state does not 
depend only on its present state. Rather, it depends on 
the past history of detections through the local oscilla- 
tor amplitude p.{t). As noted in the introduction, previ- 
ous investigations of robust unravelings in other systems 
HU have been restricted to unravelings that do have 
these properties. It is therefore interesting to ask for the 
present system, how close to the MRU is the most robust 
unraveling that is continuous and Markovian? 

Continuous Markovian unraveling (CMU) of the 
RF master equation ( |3.l| ) can be represented by a 
two-parameter family of nonlinear SSEs for the non- 
normalized state vector \il>(t)) of the form 

d\4>(t))=dt ~iH ~ Tper + J(i)er |$(t)), (6.1) 

which is to be interpreted in the ltd sense @. Here J(i) 
is a complex "current" given by 

J(t)dt = j(va + a J ')dt+ yfi dW(t), (6.2) 

where v is a complex number satisfying 



v*v < 1, 



(6.3) 



and the angle brackets denote a quantum average using 
the normalized state vector \ip(t)). We use rather 
than \tp) because the norm of \tp) has no interpretation 
in terms of probability, unlike that of \ip) in Sec. IV. The 
stochastic term dW(t) is a complex Gaussian white noise 
term satisfying 



E[dW] = 0, 
E[dW*dW] = dt, 
E[{dW) 2 ] = vdt. 



(6.4) 
(6.5) 
(6.6) 



The complex parameter v comprises the two parame- 
ters for the fam ily of unrave ling s of the form of Eq. (3.1). 
Fro m E q. (6.4) and Eq. ( |6.5D it can be shown that 
Eq. (2^6) is satisfied for all v. Thus 



p(t) = E 



m)\m) 



(6.7) 



obeys the RF master equation (3.1), while in an individ- 
ual trajectory the state remains pure. In this case there is 
no simple way to find the steady state ensemble of pure 
states. A numerical simulation of Eq. (|6.l| ), with time 
averages replacing ensemble averages, is the only way to 
proceed. 

Finding the maximally robust continuous Markovian 
unraveling (MRCMU) in this case reduces to a search 



over the ball \v\ 2 < 1 in the complex plane. We find that 
the MRCMU is for v = 1. For this value of v the "cur- 
rent" J(t) (which is real in this case) has a deterministic 
part equal to 



E [J(t)} = 7 (O 



(6.8) 



That is, the measurement yields information about the 
a x quadrature of the atomic dipole. As a consequence 
it tends to localize the atom near the a x eigenstates. 
This localization is relatively stable for high driving, since 
these are also eigenstates of the Hamiltonian fta x /2. This 
is shown in Fig. 0(e), which is a stochastically generated 
sample of 10000 states from the equilibrium ensemble for 
the MRCMU. In complete contrast to the direct detec- 
tion ensemble in Fig. 0(a), most of the states lie close 
to the a x eigenstates. For fi large as in this figure, this 
means close to the two members of the MRU shown in 
Fig. 0(c). 

Because a typical member of the MRCMU ensemble is 
fairly close to the a x axis, it is relatively little affected by 
the Hamiltonian, which causes rota tion around that axis. 
Under the RF master equation (3J) its evolution consists 
of decay to the equilibrium, with relatively small Rabi 
oscillations superimposed. This is as shown in Fig. 0(f). 
As a consequence, although the survival probability oscil- 
lates, it remains above its equilibrium value and decays 
towards it at a rate proportional to 7. This is shown 
in Fig. ^| for il = IO7. Also shown, for comparison, is 
the survival probability for the minimally robust CMU, 
which occurs for v = — 1. This is almost identical to the 
corresponding curve for direct detection in Fig. 0, as the 
ensemble is confined to the u = great circle in both 
cases. 

The survival time for the MRCMU is shown in Fig, . 
Because of our definition of survival time in Eq. ( 2.12 ) , 
the survival time for the MRCMU is determined by the 
rapid (for O 3> 7) oscillations in the survival probability 
rather than the slow mean decay. Thus it is qualitatively 
similar to the survival time for the direct detection en- 
semble, starting at 21n27~ 1 and then falling like as 
£1 increases. 



VII. COMPARISON OF THE MRU AND THE 
MRCMU 

We have seen that no continuous Markovian unraveling 
(CMU) is as robust as the maximally robust unraveling 
(MRU), which is neither continuous nor Markovian. Fur- 
thermore, the robustness for the maximally robust CMU, 
as measured by the survival time, scales in the same way 
as that for the very non-robust direct detection scheme. 
However, because of the arbitrariness in any definition of 
survival time, we are interested more in the robust en- 
sembles themselves than in the numerical values of their 
survival times. As discussed above, the MRCMU realizes 
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an ensemble that has a common feature with that real- 
ized by the MRU: the states in the ensemble tend to have 
well-defined values of u — (cr x }. In this section we wish 
to answer the question: just how close is the MRCMU 
ensemble to the MRU ensemble? 



A. Closeness of Two Ensembles 

To answer this question we require a measure (not nec- 
essarily transitive) from ensemble E A to ensemble E B , 
where each ensemble represents the same state matrix. It 
seems best to choose an operationally defined measure, 
which we do as follows. 

By allowing the same projector to reappear with dif- 
ferent indices, we can write, to any desired degree of ac- 
curacy, 



E A = {{\<j) k ),N- 1 ):l<k<N} 1 

-In 



E B = {(\^),N^):1 



< a < 



N} 



(7.1) 
(7.2) 



where \(f>k) and are normalized states such that 



N- 



(7.3) 



where N is an arbitrarily large integer. 

To define the closeness of the ensembles, imagine that 
there are two people, Alice and Bob. Alice has in her 
possession a measuring device with N settings, corre- 
sponding to the N projectors {\<fik)(4>k\ ■ k}. If setting 
k is chosen then the device makes a projective measure- 
ment with projector \4>k)(4>k\- Bob has in his possession 
an ensemble of quantum states {|^) : p}. It is Bob's 
aim to try to convince Alice that he is actually in pos- 
session of the ensemble {\(t>k) ■ k}. He must submit each 
of his N systems \ip^) to Alice, telling her which of her 
states \(j>k) each is supposed to be in. She then makes the 
appropriate measurement and, unless Bob's ensemble re- 
ally is the same as Alice's, is likely to find errors some 
of the time. An error is when a state that Bob claims is 
\4>k) is found to give the answer "no" to Alice's projective 
measurement "is the state |0fe)?" Assuming Bob chooses 
a good strategy, then the higher the probability of error, 
the larger the distance between the two ensembles. 

We can formalize this as follows. Say Bob actually 
sends state but claims it is |<^fc( M )), where the func- 
tional dependence here indicates that Bob makes his 
choice of index k based on his actual state. Then the 
probability of error for this state is 



1 



(7.4) 



The ensemble average probability of error, for Bob's op- 
timum strategy, is 



topt 



N 1 min 

M 



)Im> 



(7.5) 



where the minimum is over all one-to-one mappings M 



M 



Mm) 



M~ 



(7.6) 



Bob's strategies have to correspond to a mapping of this 
form because unless Bob names each of Alice's states 
once and once only Alice would know that Bob is lying 
when he claims to be in possession of Alice's ensemble 
{\4>k) -k). 

We could take the distance between the ensembles to 
be equal to this minimum average error probability. How- 
ever, if the state matrix p is close to being pure (with 
Tr[p 2 ] close to one), then the average error probability 
would be small regardless of what strategy Bob chose. 
In particular, if Bob had a totally random strategy then 
the error probability for an individual state would 
be, on average, 



and the ensemble average error probability would be 

(l-Tr[p|Vv}(^|]) = l-Tr[p 2 ]. (7.8) 



N 



N 

X E 



That is, the average error probability would be close to 
zero even though Bob does not take into account the dif- 
ference between his states so that his effective ensemble 
consists of N copies of the mixed state p. 

For this reason, it seems better to define the distance 
between the ensembles by the normalized error probabil- 
ity: 



d{E A \E l 



1 - min M Y^=i ~k | (gW) I VV) I 
1 - Tr[p 2 } 



(7.9) 



In this case it is easy (at least for a two-level system) to 
see that 



< d{E B \E A ) < 1, 



(7.10) 



where the lower bound is attained if and only if the en- 
sembles are identical, and where no tighter upper bound 
can be found for a given p. Thus the two ensembles could 
be said to be close if and only if d(E B \E A ) <C 1. 



B. Closeness of the MRU and MRCMU Ensembles 

In the case at hand the reference ensemble E A is the 
most robust ensemble of Sec. [v|, while the ensemble E B 
whose closeness we wish to gauge is that of the most ro- 
bust continuous Markovian unraveling of Sec. VI. The 
fact that E A has only two elements, whereas E has an 
infinitude of elements causes no problems. We simply al- 



low an arbitrarily large number TV of elements for each 
ensemble and let half of those for ensemble A be the state 
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with pro jector P + and half the state P_, as defined in 
Eq. Q. 

Because the two states P± differ only by the sign of 
(<x x ), and because both ensembles are symmetric about 
a reflection in the y — z plane, Bob's best strategy is easy 
to find. For each of his states \ipfj) he tells Alice that it 
is state P+ if it has a positive mean a x , and P- if it has 
a negative mean a x . If |VVi) is the state 

IVvKVvl = T^ 1 + u n a * + v P a v + w n a z), (7-11) 
then the probability of error for this state is 



= 1 - (^l-PsignKjIVv) 



(7.12) 



1 1 

2 ~ 2 



(M V 1 - Vss - + v fiVas + WnZ m ). (7.13) 



Averaging over all of Bob's states we get 



^SS "I" Z SB ) 1 



(7.14) 



and the distance from the MRU ensemble to the MR- 
CMU ensemble is 



d{E A \E B ) = 1 



Eflul 



Vss ~^ss 



(7.15) 



Thus to find how close the MRCMU ensemble is to 
the MRU ensemble we simply need to evaluate E[|it|], 
the ensemble average of | (a x ) | for the former. The re- 
sult is plotted in Fig. |[ The distance is always less than 
about 0.3, the value to which it appears to asymptote for 
large fl. Since this is moderately small compared to one, 
we can say that the two ensembles are moderately close. 
This is in stark contrast to either the direct detection 
ensemble or the v = — 1 CMU ensemble. For these two 
ensembles, u = for all members, so the distance to the 
MRU ensemble is the maximum value of unity. 

Also plotted in Fig. ^| is the distance from the CMU 
with v = to the MRU ensemble. This CMU has a 
number of special properties and is sometimes known as 
quantum state diffusion (2^j2^] . We see that it also gives 
an ensemble whose distance to the MRU ensemble is less 
than unity. However, the distance is considerably greater 
than that of the MRCMU, asymptoting to a value greater 
than 0.4. 



VIII. DISCUSSION 

The resonance fluorescence master equation generates 
surprisingly rich dynamics for a two-level atom. Here we 
have investigated those dynamics using the technique of 
finding the maximally robust unraveling. That is, finding 
the scheme for monitoring the fluorescent radiation that 



collapses the atom into pure states that are, on average, 
the most robust. By robust we mean that they survive 
best under the master equation evolution so that, once 
the monitoring ceases, the probability for the atom to 
be found at some later time to still be in the state into 
which it was collapsed by the monitoring, is maximized. 

The property of producing robust states may give 
the maximally robust unraveling potential applications 
particularly in quantum information technology p3|, 24 
where minimizing the effect of environmental decohcr 
ence is essential. Quite separately from any application, 
the maximally robust unraveling is useful for character- 
izing the dynamics of open quantum systems, such as 
the resonantly driven atom of this study. It has been 
suggested before J|,[| that the ensemble arising from the 
MRU is the most natural representation of the system's 
stationary state matrix in terms of pure states (state vec- 
tors). 

In this work we have found that the MRU for the RF 
master equation is an adaptive interferometric monitor- 
ing scheme proposed in Ref. p0| . The atom's radiation 
is, prior to detection by a photodetector, interfered at a 
beam splitter with a reflected local oscillator. The mea- 
surement is adaptive because the local oscillator ampli- 
tude (which is comparable in magnitude to the field ra- 
diated by the atom) has its phase changed by it every 
time a detection occurs. This detection scheme has the 
remarkable property that, in steady state, the atom sim- 
ply jumps between two fixed pure states. In the large 
driving limit, these two pure states are close to eigen- 
states of the driving Hamiltonian Qa x /2. 

The adaptive interferometric monitoring scheme was 
designed in Ref. |^0| specifically to have this property 
of producing a stationary ensemble containing just two 
members. In that reference it was found that other de- 
tection schemes, such as spectral detection (resolving the 
three Mollow peaks), and another adaptive scheme, give 
rise to similar behaviour. That is, the atom jumped be- 
tween states that were close to a x eigenstates in the large 
driving limit. In Ref. (2(]] it was speculated that this be- 
haviour was what the atom "wanted to do". Here we 
have confirmed that the resulting two-member ensemble 
is indeed the most robust and hence arguably the most 
natural. In the context of the study of decoherence and 
the classical limit, it appears that jumping between two 
fixed states is the most classical behaviour for a strongly- 
driven two level atom. 

Another issue we have investigated in this paper is how 
close to the MRU one can approach if one restricts the 
unravelings to continuous Markovian ones. In the con- 
text of the fluorescent atom, this means unravelings re- 
alizable from homodyne measurements. They give rise 
to evolution on the Bloch sphere that is continuous (but 
not differentiable) and Markovian. This is an interesting 
question because the set of continuous Markovian unrav- 
elings is easily parameterized by real numbers, unlike the 
set of all possible unravelings, which is too large to be 
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finitely parameterized in this way. For this reason, pre- 
vious work in MRU J||| has concentrated on finding the 
most robust CMU. 

In this work we have found that the MRCMU has a 
robustness, as measured by the survival time, which falls 
as fl^ 1 as the driving increases. This is similar to the 
result for direct detection, and contrary to that of the 
MRU for which the survival time is constant at 2 In 27 _1 . 
However, as we have shown graphically, the distribution 
of states on the Bloch sphere for the MRCMU ensem- 
ble is qualitatively much closer to that of the MRU than 
to that of direct detection. Furthermore, we have in- 
troduced a quantitative measure for the closeness of two 
ensembles of pure states, and applied this to the vari- 
ous unravelings. We find that the MRCMU ensemble is 
reasonably close to the MRU ensemble, while the direct 
detection ensemble is as distant as is possible from the 
most robust ensemble. 

This result has wider implications in the program of 
decoherence, robustness, and the classical limit. A two- 
level atom is an extremely non-classical system. In the 
limit of strong driving the stationary state matrix is al- 
most fully mixed, and the existence of two discrete levels 
manifests strongly in the dynamics: the maximally ro- 
bust unraveling has jumps between two almost orthogo- 
nal states. By contrast, under a continuous Markovian 
unraveling the atomic state does not jump, but rather 
diffuses around the Bloch sphere. 

We have shown that, despite the atom's nonclassical- 
ity the ensemble generated by a CMU can be reasonably 
close to that generated by the MRU. This suggests that 
restricting an investigation to continuous Markovian un- 
ravelings is not a serious restriction (provided that one 
is not interested in the numerical value of the maximum 
survival time) . The basic nature of the maximally robust 
unraveling should be evident from that of the maximally 
robust CMU. For the two-level atom it is fortunate that 
the absolute maximally robust unraveling can be found 
analytically. For more general systems a numerical search 
would be necessary, and, to be practical, would have to 
be confined to finitely parameterizable unravelings such 
as the continuous Markovian unravelings. Thus our re- 
sults lends credence to the whole program of finding the 
(approximately) maximally robust unraveling for open 
quantum systems. 



Now what we want is p(t), the probability that the last 
detection was at a time t before the present, at steady 
state. That is, the probability that there was a detection 
at time t in the past, and that there were no detections 
from then until now. In other words, 



p(t) = P[n(t) A d(t)}, 



(A2) 



where P[A A B] means the joint probability of A and B. 

Now from the definition of conditional probability, 
P[AA B] = P[A\B]P[B]. Therefore 



p(t) = P (t)P[d(t)}. 



(A3) 



But at steady state (that is, after initial transients have 
decayed), the probability that there was a detection at 
time t in the past, given no other information, does not 
depend on t. That is, P[d(t)] is a constant, so we simply 
have 



P (t)=M~ 1 P (t), 



(A4) 



for some constant N '. 

To find the constant of proportionality we just note 
that, since the last detection must have been at some 
time t in the past, J °° p{t)dt = 1. Here we are actually 
treating p{t)dt as the probability that the last detection 
was in the interval [t, t + dt) . From this condition it is 
easy to see that 



/>oo 

N= / P (t)dt, 

Jo 



(A5) 



which gives Eq. (4.7) 
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APPENDIX: DERIVATION OF EQ. (4.7) 

Let n(t) be the event that there were no detections 
from time to to the present time to + t. Let d(t) be the 
event that there was a detection at time t before the 
present. Then 



P (t) = P[n(t)\d(t)], 



(Al) 



where Po(t) is as in Eq. (4.6) and where P[A|_B] means 
the probability of A given B. 
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Figure 1 is attached at the end. 

FIG. 1. Plots (a), (c), and (e) show the distribution of pure 
states in the stationary ensemble on the Bloch sphere under 
various detection schemes. Plots (b), (d), and (f) show the 
decay towards the stationary state matrix of a typical member 
of each ensemble. Plots (a) and (b) are for direct detection, 
(c) and (d) for the maximally robust unraveling, and (e) and 
(f) for the maximally robust continuous Markovian unravel- 
ing. For all plots fi/7 = 10. The ensemble in (c) consists 
of just two members, and the dots indicating their positions 
have been enlarged to make them more easily visible. 
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Time 



FIG. 2. Plot of the ensemble average survival probability 
S(t) versus time (in units of 7 -1 ) for the stationary ensemble 
of direct detection. The three curves are, from the top down, 
n = 7/2, n = 7, and n = IO7. 



FIG. 3. Plot of the ensemble survival time (in units of 7 -1 ) 
versus fi/7. The solid curve is for the ensemble arising from 
direct detection (Sec. IV). The dashe d cur ve is the large fi/7 
analytical approximation to it in Eq. (4.15). The dotted curve 
is that from the maximally robust unraveling (Sec. V). The 
dash-dot curve is that from the maximally robust continuous 
Markovian unraveling (Sec. VI). The unevenness in this final 
curve is due to statistical error in the ensemble averages V v 
and V w ), and gives an indication of the size of the statistical 
error. 
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FIG. 4. Diagram of the experimental configuration for an 
adaptive interferometric measurement of the fluorescence of 
an atom. The signal and a local oscillator are incident on a 
beam splitter of reflectance r<l. The amplitude of the local 
oscillator is variable as a function of time, determined by an 
electro-optic modulator (EOM). The modulator is controlled 
by the experimenter, being inverted every time a detection 
occurs. 
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FIG. 5. Plot of the ensemble average survival probability 
S(t) versus time (in units of 7~ ) for the stationary ensem- 
bles of various detection schemes all with fi/7 = 10. The 
three curves are for the minimally robust (v = — 1) continuous 
Markovian unraveling (solid), the maximally robust (v — 1) 
continuous Markovian unraveling (dash-dot) and the maxi- 
mally robust unraveling (dotted). 
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FIG. 6. Plot of the distance of various ensembles from the 
maximally robust ensemble as a function of fi/7. The solid 
line is for direct detection, which is the same as for continu- 
ous Markovian unraveling (CMU) with v = — 1. The dashed 
line is for a CMU with v = 0. The dash-dot line is for v = 1 
(the maximally robust CMU). In the latter two cases, the 
three lines indicate the mean, and the mean plus or minus 
one standard deviation (the statistical error in the ensemble 
average) . 
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